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We investigate the impact of the neutron-skin thickness, AR„p, on the energy difference between 
the anti-analog giant dipole resonance (AGDR), -Eagdr, and the isobaric analog state (IAS), -Bias, 
in a heavy nucleus such as ^°®Pb. For guidance, we first develop a simple and analytic, yet physical, 
approach based on the Droplet Model that linearly connects the energy difference Bagdr —Bias with 
ARnp. To test this correlation on more fundamental grounds, we employ a family of systematically 
varied Skyrme energy density functionals where variations on the value of the symmetry energy 
at saturation density J are explored. The calculations have been performed within the fully self- 
consistent Hartree-Fock (HF) plus charge-exchange random phase approximation (RPA) framework. 

We confirm the linear correlation within our microscopic apporach and, by comparing our results 
with available experimental data in ^°®Pb, we find that our analysis is consistent with ARp„ = 0.204 
± 0.009 fm, J = 31.4 ± 0.5 MeV and a slope parameter of the symmetry energy at saturation of 
L = 76.4 ± 5.4 MeV — the attached errors correspond to a lower-limit estimate of the systematic 
plus experimental uncertainties. These results are in agreement with those extracted from different 
experimental data albeit, L and ARpn, are somewhat large when compared to previous estimations 
based on giant resonance studies. 
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I. INTRODUCTION 

Different experimental methods, either direct or in¬ 
direct, have been proposed to extract the value of the 
neutron-skin thickness in finite nuclei, that is, the dif¬ 
ference between neutron and proton root-mean-square 
radii, 

^ ( 1 ) 

The neutron skin thickness is an observable that has 
kept much attention from both experimental and the¬ 
oretical viewpoints. This is because it is one of the 
most promising observables in nuclear structure to con¬ 
strain the density dependence of the symmetry energy 
around the nuclear saturation density [IM3- The sym¬ 
metry energy plays an important role in understand¬ 
ing the mechanisms of different phenomena in nuclear 
physics and nuclear astrophysics [sl-l^: it directly af¬ 
fects the properties of exotic nuclei, the dynamics of 
heavy-ion collisions, the structure of neutron stars, and 
the simulations of core-collapse supernova. 

The Lead Radius Experiment (PREX) at the 
Jefferson Laboratory has provided the first model- 
independent evidence on the existence of a neutron- 
rich skin in ^°®Pb [2^. Relying on the fact that the 
weak charge of the neutron is much larger than the 
corresponding proton one, PREX used parity-violating 
electron scattering to probe the neutron distribution of 
^®®Pb. To foster this field, more experiments have been 


already approved with both, the aim of improving the 
reached accuracy in ^°®Pb and to explore other mass 
regions. On the other side, neutron densities have been 
traditionally probed mostly by nucleon or a scattering. 
For example, by using proton elastic scattering on Sn 
and Pb isotopes [s^; or by measuring photons emitted 
during the decay of antiproton states [3l|, • One can 

also obtain information on the neutron skin thickness 
from giant resonance properties, such as the excitation 
energy of the isovector giant dipole resonance (IVGDR), 
the total electric dipole polarizability (q:d), the excita¬ 
tion energy of the isovector giant quadrupole resonance 
(IVGQR) or, yet with more warnings, from the energy 
and strength of the pygmy dipole resonance (PDR) in 
neutron-rich nuclei Last but not least, the total 

strength of the charge-exchange spin-dipole resonances 
(SDR) can be related to the neutron skin in a very trans¬ 
parent way [ 3 , . It is important to mention, how¬ 

ever, that all hadronic probes require model assump¬ 
tions to deal with the strong force introducing possible 
systematic uncertainties. 


Recently, the authors of Refs, have proposed 

a new method to extract the neutron skin thickness 
based on the measurement of the excitation energy of 
the anti-analog giant dipole resonance (AGDR), that 
can be observed in the charge-exchange (p, n) reaction. 
The AGDR was first studied experimentally in Ref. [55| . 
Already in Ref. [s^, the authors had pointed out that 
the excitation energy of the AGDR is sensitive to the 
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neutron-skin thickness. More recently, the energy differ¬ 
ence between the AGDR and the isobaric analog state 
(IAS), i^AGDR — Aias, in ^'’®Pb has been obtaine d by 
measuring the direct 7 -decay between these states . 

In this paper, we shall analyze the relationship of 
the neutron-skin thickness and the energy difference 
Aagdr — Bias , by using a fully self-consistent Hartree- 
Fock (HF) plus charge-exchange random phase approx¬ 
imation (RPA) framework with a family of Skyrme en¬ 
ergy density functionals. We try to understand also the 
qualitative features of such relationship through a sim¬ 
ple, yet physical and transparent model. Our approach, 
as compared to Ref. , incorporates specific effects of 
the Skyrme functionals such as the effective mass and 
the isovector enhancement factor (cf. Sec. IIIII) . Then, 
by comparing the theoretical and experimental results 
for Aagdr ~ Bias-, we extract the neutron-skin thick¬ 
ness in ^°®Pb. This allow us, in turn, to estimate the 
compatible values for the symmetry energy J and its 
slope parameter L (at nuclear matter saturation den¬ 
sity). The extracted values of J and L are eventually 
compared to the results obtained by other analysis on 
different observables. 

The outline of the paper is the following. In Sec. |lT] 
the theoretical model is briefly presented: we focus, in 
particularly on the charge-exchange random phase ap¬ 
proximation (RPA) based on the use of non-relativistic 
Skyrme energy density functional (EDFs). In Sec. IIIII 
we derive our analytic model to explain the relationship 
between the energy difference Aagdr — Bias of AGDR 
and IAS, and the neutron-skin thickness. A detailed 
quantitative analysis of such correlation is performed 
by employing a family of so-called SAMi-J Skyrme func¬ 
tionals, in Sec. 113 In Sec. El we summarize the results 
and draw our conclusions. 


II. MICROSCOPIC MODEL: 

CHARGE-EXCHANGE RPA 

The calculations are done within the framework of the 
Skyrme HF plus charge-exchange RPA. We adopt 
the standard form of Skyrme interactions with the no¬ 
tations of Ref. [ 5 ^ . Two nucleons characterized by the 
space, spin and isospin variables r^, cr^ and Ti interact 
through a zero-range, velocity-dependent and density- 
dependent force that reads 


H(ri,r-2) 

= 

to(l + xoPa)5[r) 


-f 

]^ti{l+xiPa)[P'^5{r) + 5{r)P‘^] 


-h 

^ 2(1 + X2Pa)P' ■ 5{r)P 


-f 

+ X3Pa)p'^{R)5{r) 

0 


-h 

iWo{ai -1- 0 - 2 ) • [P' X (5(r)P] , 1 

where r = ri 

-r2, R= \{ri + 72), P = ^(Vi - V 


P' is the hermitian conjugate of P (acting on the left), 
Pa = + (Ti ■ < 72 ) is the spin-exchange operator, and 

P = Pn+Pp is the total nucleon density. Within the stan¬ 
dard formalism, the total binding energy of a nucleus 


can be expressed as the integral of the Skyrme den¬ 
sity functional which includes the kinetic-energy 
term /C, a zero-range term T-Lq, the density-dependent 
term "Ha, an effective-mass term Tleff, a- momentum 
dependent term (that mimics finite-range effects) 'Hfin, 
a spin-orbit term T-Lso, a spin-gradient term T-Lsg, and a 
Coulomb term Hcoui- 

Here, we will briefly summarize the formulas for the 
charge-exchange RPA calculations. The well-known 
RPA method |CT . | 6 l| in matrix form is given by 


f A B 
\B* A* 
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where B^, is the energy of the i^-th charge-exchange 
RPA state and are the corresponding forward 

and backward amplitudes, respectively. The matrix el¬ 
ements A and B are expressed as 

Ami,nj — (^m T (mjj I^-es j (4) 


Bmi,nj — (^^|Pres|G)- (3) 

where the sub-indexes i,j refer to occupied states, m,n 
to unoccupied satates, VJ-es is the residual interaction 
and e the single-particle states energy. 

The particle-hole (p-h) matrix elements are obtained 
from the Skyrme energy density functional including all 
the terms (the Coulomb term T-Lcoui is not active in this 
case). The explicit forms of the matrices A and B are 
given in Ref. [62j | in the case of a Skyrme force. 

We will use the following operator for the AGDR ex¬ 
citation. 


d± = '^rABim{ri)t± , ( 6 ) 

i 

which corresponds to the AJ = 1, AL = I, AS" = 0, 
J'^ = 1“ resonance. We will also calculate the IAS. The 
IAS excitation operator reads 

OiAS = (7) 

i 

and corresponds the A J = 0, AL = 0, AS = 0, = 1+ 

states. 


In Fig. 1, we illustrate the relevant states of a tar¬ 
get and its daughter nucleus, namely the ground state, 
the giant dipole resonance and the giant Ml state of a 
target nucleus. Tq 0 labels the ground-state isospin 
of the target nucleus, that is, {N — Z)I2. The cor¬ 
responding resonance states in the daughter nucleus 
reached by the (p, n) charge-exchange reaction are also 
displayed, namely the IAS (isospin = Tq), and the anti¬ 
analog (isospin = To — 1) states: Gamow-Teller reso¬ 
nance (GTR), SDR, and AGDR. As shown in the fig¬ 
ure, the AGDR corresponds to AJ'^ = I“, AL = 1, 
and AS" = 0 excitation, and represents the anti-analog 
giant dipole resonance because it is the Tq — 1 compo¬ 
nent of the charge-exchange of the isovector giant dipole 
resonance (IVGDR). 
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constant k for a dipole excitation is determined by the 
self-consistent condition between the vibrating potential 
and density: 
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T—0 

A =1 
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( 11 ) 

( 12 ) 


for an isoscalar (r = 0) or isovector (r = 1) dipole (A = 
1 ) excitations, respectively, hioo is the major shell gap, 
« 41 MeV, Vi is the strength of Lane potential 

U = Vi and (r^) is the mean square radius. The 
classical EWSR of the isovector dipole mode is given by 


FIG. 1: Various states related to the target and daughter 
nucleus, including the ground state, the IVGDR and Ml 
excitations of the target nucleus (isospin = To), and the IAS 
(isospin = To), anti-analog states (isospin = To — 1) in the 
daughter nucleus excited in a (p,n) reaction. 


III. GUIDELINE FROM A SIMPLE ANALYTIC 
MODEL 


m„(r = l,A=l) = TAyl. (13) 

It has to be noted that the energy of the isoscalar 

dipole mode {viz., the spurious center-of-mass mode) 
turns out to be at zero energy, as it should be. This 
can be checked from Eqs. m and m- Moreover, the 
energy of the isovector giant dipole resonance, = 

Aivgdr, becomes 


In this Section we will develop a simple, yet physically 
sound, model for the excitation energies of the IVGDR, 
the IAS, and the AGDR. This effort will allow us to 
gain a deeper insight into the relevant “macroscopic” 
physics of our problem, namely the relationship between 
the neutron skin thickness and the energy difference of 
the AGDR and the IAS. 

We start from the RPA dispersion relation for a sep¬ 
arable interaction of the type kO^O, that is, 

\{p\d\h)\^ ^ \{p\d\h)\^ 1 

^ Eph- E eph + E K 


In this equation, Sph is the unperturbed p-h excita¬ 
tion energy; if we assume that there is only one unper¬ 
turbed configuration at energy Sph = £o that exhausts 
the whole unperturbed strength Sq, then the equation 
becomes 




-f 


So 


£0 — E £o -f- E 


1 

5 

K 


and, therefore. 


E^ = £q + 2ktoo 


(9) 


( 10 ) 


where mg = EqSq is the energy weighted sum rule 
(EWSR). 

For each mode, there will be a different coupling con¬ 
stant; however, the isospin invariance dictates that the 
coupling constants of the AGDR and the non-charge ex¬ 
change IVGDR to be the same. We start by considering 
the well-known IVGDR case: O = DEio(u)Ti*^ — 
note that for the macroscopic model we use the compo¬ 
nents of r instead of t. 

We will use the Bohr-Mottelson quantal harmonic os¬ 
cillator (QHO) model [i^; in this case, the coupling 


^^IVGDR - So (l + > (14) 

and it is well known that by using standard values for 
the unperturbed energy (eq = fuvo = 4I/A^/^ MeV) as 
well as for the radius ((r^) = 3Rq/5 = 3(1.2)^A^/^/5 
fm^), together with Vi = 130 MeV, the excitation en¬ 
ergy provided by Eq. (ITTl) is EJzl 80/A^/^ MeV, in 
good agreement with the empirical systematics for the 
IVGDR in the mass region A > 40. 

There are two important differences between the in¬ 
teraction assumed in the Bohr-Mottelson model and the 
Skyrme interaction. In the former case, the effective 
mass is taken to be m*/m = 1, while for Skyrme in¬ 
teractions this value depends on the chosen set, being 
in uniform matter as well as in the interior of nuclei 
close to the empirical value, m* jm w 0.7. The effective 
mass changes the unperturbed energy to be e = eo 
Moreover, in the case of momentum-dependent inter¬ 
actions such as the Skyrme forces (or other non-local 
forces), the classical EWSR should be multiplied by 
1 + a where a is the so-called enhancement factor. For 
the dipole case, a is typically around ~ 0.2 [o^. Gon- 
sequently, Eq. da should turn into 

E"^ = -^Eg + 2K.'mg (1 -I- a). (15) 

It is a simple exercise to show that Eq. m reproduces 
the experimental IVGDR systematics as well as Eq. da 
if the coupling constant k' is reduced with respect to k, 
that is, k' ~ 0.7k. By inspecting Eq. da we can 
also conclude that this implies a quenched value for the 
strength of the Lane potential V{ with V{ w 0.7Ul (this 
value is of course indicative, in keeping with the rather 
crude approximations of this analytic model). 

In the case of charge-exchange excitations of nuclei 
having a neutron excess, like ^°®Pb, the Tamm-Dancoff 
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TABLE I: Excitation energies and EWSRs (mo) of the AGDR and IVGDR for the family of SAMi-J interactions are 
displayed, together with the excitation energy of the IAS, the enhancement factor 7 = /?/(! + «) and the prediction 
of the macroscopic model -Eagdr estimated from Eq. (1^ . To evaluate -Eagdrj we take e = 41A \Jmlm* MeV, 
AEc = 2 (3/5)®^^ and the energies of the current table. See the text for more details. 


Force 

Aagdr 

[MeV] 

[fm^l 

l^IYGDR 

[MeV] 

[fm^] 

Aias 

[MeV] 

7 

7 -tmac 

-^AGDR 

[MeV] 

SAMi-J27 

28.20 

11153.99 

13.99 

4073.4 

18.89 

0.369 

25.51 

SAMi-J28 

27.94 

11119.60 

13.91 

4071.0 

18.74 

0.366 

25.25 

SAMi-J29 

27.64 

11079.18 

13.74 

4067.9 

18.58 

0.362 

24.87 

SAMi-J30 

27.37 

11055.00 

13.58 

4065.4 

18.42 

0.360 

24.51 

SAMi-J31 

27.16 

11047.33 

13.42 

4063.5 

18.28 

0.359 

24.18 

SAMi-J32 

26.94 

11055.63 

13.27 

4062.6 

18.16 

0.361 

23.89 

SAMi-J33 

26.73 

11078.24 

13.13 

4061.6 

18.05 

0.364 

23.62 

SAMi-J34 

26.54 

11113.27 

13.00 

4061.3 

17.96 

0.368 

23.40 

SAMi-J35 

26.39 

11158.82 

12.88 

4061.4 

17.89 

0.374 

23.21 


approximation (TDA) is known to provide results that 
are quite similar to those from the RPA. This is because 
the coupling between the r_ excitations and r_|_ excita¬ 
tions is small due to their quite large energy difference. 
We simplify thus the RPA dispersion relation to be the 
TDA one, 

T,h ^ ^ 


The AGDR has AL = 1, /S.S = 0, where the corre¬ 
sponding operator is 

'^riYioT^{i). ( 21 ) 

i 

The NEWSR reads 

5o(r-, A = 1) - 5o(r+, A = 1) = (22) 


As we have done already, we assume that there is only 
one unperturbed configuration at energy Sph = Sq ex¬ 
hausting the whole unperturbed strength So, so that 


gp 1 

Eq — E K 


(17) 


The solution of this simplified TDA equation is 


E 


Eq -h kSq 
Too 


Eq T 


£0 


(18) 


Let us consider the IAS first. It is well known that, to 
a first approximation, its excitation energy is associated 
with the Coulomb energy shift AEc between the parent 
and the daugther nuclei. In our TDA model we can 
write the unperturbed p — h energy as eq = —U + AEc 
where U is, as above, the Lane potential. The non¬ 
energy weighted sum rule (NEWSR) obtained by using 
the operator r_(i) is 2(iV — Z). Therefore, 

Eias = -U + AEc + Kl=j2(iV - Z) (19) 


and if 



( 20 ) 


the IAS energy coincides with AEc- 

Let us finally move to the subject of main interest 
for us, namely the AGDR. Our goal is to have a trans¬ 
parent interpretation of the results obtained with the 
microscopic Skyrme model. 


where 


_ - Z{t'^)p 

N - Z 


(23) 


In this sum rule the t_ contribution is largely dominant 
in nuclei with neutron excess like ^°®Pb; the same dom¬ 
inance holds for the energy-weighted sum rule, that can 
be written as 


3 h^A 

mo(r_, A = 1) - mo(r+, A = 1) = —- — {l + a + P). 

ZTT 2m 

(24) 

a is the same as in the IVGDR case that we have dis¬ 
cussed above, whereas the definition of j3 can be found, 
in the case of a Skyrme interaction, in Ref. 

Within the framework of our approximation, the 
AGDR unperturbed energy can be written as e — [/ -I- 
AEc', consequently, its TDA energy from the simplified 
equation (fTSl) reads 


l^AGDR 


E-U + AEc + 


V{ {N{r^)n - Z{r'^)p) 

2 A(r 2 ) 


E-U + AEc + 


VI 1^(1 + « + ^) 

{r^) E — U + AEc 


(25) 


We are supposed to use the same coupling constant that 
has been already used in the case of the IVGDR; ac¬ 
cording to our previous discussion, this will be different 
from the case of the Bohr-Mottelson model if used in 
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conjunction either with an effective mass and/or with 
an enhancement factor as in the Skyrme case. For con¬ 
venience, we shall define here Vi = 1^/(1 + a) — note 
that Vi = Vi if m* jm = 1 and V\ ~ Vi even for realis¬ 
tic models with m* /m < 1. One can also notice that a 
simplification of Eq. (l25|) comes from the fact that for a 
heavy nucleus such as replacing e—U+AEc with 

AEc will produce an error of only a few %. Specifically, 
if we assume m* jm, w 0.7 and V\ « 130 MeV as previ¬ 
ously done, e — U = ^mjm* — V[(N — Z) j^A Ri 

1.3 MeV which correspond to about 7% when compared 
to AEc- We use this simplification to write the energy 
difference between the AGDR and the IAS: 


Aagdr — Aias 


Vi 

2 ( 1 -fa) 


- 2 ’(r 2 )p) 

A{r‘^) 


(r^) AEc{la) 


(26) 


It is convenient to define the quantity 7 = /3/(l -f a) 
since it is almost constant if we consider the interac¬ 
tions employed in the current study (see Table |l|. Fi¬ 
nally, approximating the IAS energy as the Coulomb 
shift energy between parent and daugther nuclei by 
AEc = 2 (3/5)^^^ Z!, we may write 

Vi(i- 1 - 7)3 
AEc 2 2 mc^(r^) 

5 /5 Vi(l-f7) he 
sV 3 a^iZ 


Aagdr — Aias = 


If we take a ~ 0.2, 7 k- 0.4 and Vi ~ 130 MeV, we 
find Aagdr — Aias ~ 9 MeV, which is in reasonable 
agreement with the result of our realistic calculations. 

This schematic model gives us the opportunity to un¬ 
derstand the sensitivity of Aagdr —Ejas on the neutron 
skin thickness Ai?„p. In fact, as it was done in Ref. 
to which we confer the reader for details, we can relate 
the interaction strength of the potential Vi with the 
neutron skin thickness via the Droplet Model (DM), 


Vi^8lasym(A)-eF^/3]. (28) 


The DM also predicts that 

7 3J 1 

- 2{r^y/^ I - Ic 

X (^AR„p - ARlf + PcirY^^^ , (29) 


where Ic = e^Z/2QJR is a Coulomb correction to the 
total neutron excess / = {N — Z)IA, asym(A) is the 
symmetry energy parameter of the DM, and is 

a surface correction to the neutron skin thickness due 
to the different neutron and proton surface diffuseness. 
The latter quantity has been shown to be approximately 
constant in 2°®Pb (AR^”'" w 0.09 ± 0.01 fm) when cal¬ 
culated by a large set of energy density functionals of 
different kind | 66 | . Since Ic corresponds to a correction 


of about a 10 % to / in heavy neutron-rich nuclei such as 
^°®Pb, we will assume in what follows that I — Ic ~ I, 
and find 


J (^sym(A) 


3J ^ ARnp - AR™"f 3 Ic 

2J ^ p)i75 


(30) 


and by combining this result with Eq. (113 one finds 
that 


Eagdr — Aias 


g / 5 J 1 -f 7 he 


3 / UYiZ 




V 3J 



( AR, 


■np 


, 2 ) 1/2 I 




(31) 


For a given nucleus, Eq. (ISTl) predicts an explicit linear 
anti-correlation of Aagdr — Eias with ARnp. We will 
show in the next Section that this correlation is actually 
displayed by the microscopic results. 

We have also found very instructive to relate the dif¬ 
ferent excitation energies within our macroscopic model, 
and check if the microscopic results follow such relation¬ 
ship. In doing that we have used the TDA expressions 
for the IVCDR, IAS and ACDR, and after some alge¬ 
bra, we arrive at 


Eagdr 


AEc 



e-U\ 

A^j 


(32) 


(Eivgdr — e)2(l + 7) 


e 

AEc 1 -I- 


1 


s-U 

AEc 


As previously done, in a nucleus such as ^°®Pb, (e — 
U)/AEc can be neglected. Therefore, within a good 
approximation, we can write 

Eagdr - Eias « . ^ (Eivgdr - e) 2(1 -|- 7 ) (33) 

IXtLc 

or 


S ^AGDR 

Eagdr - Eias « (Eivgdr - e) ^ivgdr ■ (^4) 


We define the energy of the ACDR extracted from Eq. 
(IMl) as E™qqj^, in Table I and hereafter. For the SAMi- 
J family and for a fixed nucleus, this formula suggests 
that the energy difference Eagdr — Eias should display 
the same trends as shown by Eivgdr- In fact, e de¬ 
pends only on the effective mass which is constant for 
the SAMi-J family (e = A\A~^I'^\Jmlm* MeV). More¬ 
over, AEc is expected to not vary, and 1 -I- 7 is also ap¬ 
proximately constant (cf. Table |I|. The expression (|M)) 
reflects the idea that the physics encoded in the energy 
difference Eagdr — Eias reflects that of the IVCDR, as 
expected because of isospin invariance. 

In Table U we present the predictions of the SAMi- 
J family for the different observables under study. The 
reader can verify that the latter equations of this Section 
reasonably reproduce the microscopic HF-RPA results: 
although there is an almost constant shift, the trend of 
Eagdr (R-PA), however, is almost perfectly reproduced 
by the value E™qqj^ in Eq. (IMl) . This finding gives us 
confidence in using the simple arguments in this Section 
to interpret the microscopic results. 
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IV. RESULTS AND DISCUSSIONS 

In this Section, we discuss the results obtained by 
employing the SAMi-J Skyrme energy density function¬ 
als to calculate the HF ground state and RPA excited 
states. The SAMi-J interactions are characterized by 
different values of the symmetry energy at saturation 
density: this value varies between 27 MeV and 35 MeV 
(in steps of 1 MeV), and the force parameters are fitted 
using properties of selected nuclei while keeping at the 
same time the constraints on few properties of nuclear 
matter (nuclear incompressibility Koo = 245 MeV, and 
nucleon effective mass m*/m = 0.675). For details, the 
reader should consult Ref. [4l |. 

The ground state properties of ^°®Pb are calculated in 
coordinate space using box boundary conditions. The 
radius of the box is taken to be 20 fm: the same box is 
used to calculate discrete states at positive energy that 
are associated with the continuum part of the spectrum. 
A cutoff energy of 60 MeV (in the single-particle energy) 
is adopted for the RPA calculations. With this energy 
cutoff, we have checked that the non-energy weighted 
sum rules for both AGDR and IAS are satisfied at the 
level of about 99.97% for all Skyrme functionals used in 
the present study. 

In Fig. 2, we show the response functions correspond¬ 
ing to the IAS and AGDR operators obtained for ^°®Pb 
by using the SAMi-J Skyrme functionals: the RPA re¬ 
sults have been smeared out by using Lorentzian func¬ 
tions. As we can see, the IAS the peak energy has small 
fluctuations as it varies between 17.5 and 18.6 MeV 
for the different SAMi-J parameter sets. As for the 
AGDR case, the peak energies vary between 26 MeV 
and 28.5 MeV by using the different SAMi-J parameter 
sets. Experimentally, the mean AGDR energy has been 
extracted from the response function in the energy in¬ 
terval 5-15 MeV above the IAS energy. To compare our 
results with the experimental findings, we shall use the 
same energy range to calculate the mean energy from 
the AGDR response. 

In Fig. 3 we display the excitation energy of the 
AGDR and IAS as a function of the nucleon effective 
mass m* jm, calculated with the SAMi-J and SAMi-m 
Skyrme functionals. We remind that the main differ¬ 
ence between the SAMi-m and SAMi-J functionals is 
that in the former case the nucleon effective mass varies 
(in steps of 0.05) when fitting the parameters while Koo, 
J and L are kept constant (as above, we refer to [l^ for 
details). The red squares in the panels correspond to re¬ 
sults from SAMi-m. The results obtained from SAMi- 
J (black circles) are displayed in such a way that the 
bottom (top) point corresponds to the highest (lowest) 
value of J. The conclusion from these panels is that the 
excitation energy of the AGDR is sensitive to the sym¬ 
metry energy at saturation density while, as expected, 
the variation of the AGDR excitation energy within the 
sets of the family SAMi-m is small. In the case of the 
IAS, the excitation energy is neither sensitive to the 
symmetry energy nor to the effective mass. 

The calculated energy differences Aagdr — ^^lAS be- 



FIG. 2: (Color online) The (a) IAS and (b) AGDR response 
functions calculated by using the SAMi-J Skyrme energy 
density functionals. The discrete RPA peaks have been 
smeared out by using Lorentzian functions with (a) 300 keV 
and (b) 3 MeV width. 


TABLE II: The values of the neutron-skin thickness of ^'^®Pb 
obtained in the present work are compared to other values 
extracted by means of different experimental methods. 


Method 

Ref. 

Date 

ARpn{fm) 

antiproton absorption 

m 

2001 

0.180 ± 0.030 

(q, a') IVGDR 

m 

2004 

0.120 ± 0.070 

PDR 

[43] 

2010 

0.194 ± 0.024 

{P,P') 

[3^ 

2011 

0.156 ± 0.025 

Q-d 

m 

2012 

0.168 ± 0.022 

parity violation 

[29] 

2012 

0.330 ± 0.170 

(7,^°) 

m 

2014 

0.150 ± 0.030 

AGDR 

present 

2015 

0.204 ± 0.009 


tween AGDR and IAS, obtained by employing the 
SAMi-J Skyrme functionals, are displayed as a func¬ 
tion of the corresponding neutron-skin thickness in Fig. 
4: in particular, the solid circles correspond to the sets 
SAMi-J27 to SAMi-J35, from left to right. As we men¬ 
tioned above, for the excitation energy of the AGDR 
we take the centroid of the theoretical strength distri¬ 
bution, calculated in the energy interval from 5 to 15 
MeV above the IAS energy. The results show that the 
energy differences Aagdr — Eias between AGDR and 
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FIG. 3: (Color online) Excitation energies of the IAS and 
AGDR in as a function of effective mass and symmetry 
energy at saturation density. The filled boxes are the results 
of the SAMi-m family, while the filled circles are those of the 
SAMi-J family. The top circle in the line of circles of each 
window corresponds to the lowest J value (J = 27 MeV) and 
by going down one increases the J value in steps of 1 MeV. 
The lowest circle then corrsponds to the maximum value, J 
= 35 MeV. 



FIG. 4: (Color online) The energy difference Eagdr — Eias 
of AGDR and IAS as a function of neutron-skin thickness, 
obtained by using the SAMi-J family of Skyrme functionals 
consistently. The calculated values are presented as solid 
circles. Two different experimental data [HJ, [s^l are also 
shown as solid (magenta) and dashed (blue) lines, respec¬ 
tively. The arrows indicate the neutron skin constrained by 
these experimental data. 


IAS decrease with increasing values of the neutron-skin 
thickness, and a strong linear correlation exists; this is 
quite well justified by the model that has been devel¬ 
oped in Sec. uni 

In Fig. 4 we also super-impose two different experi¬ 
mental data. In Ref. (that will be denoted as Expl 
hereafter), the AGDR has been separated from other ex¬ 
citations by means of the multipole decomposition anal¬ 
ysis of the ^°®Pb(p, n) reaction at a bombarding energy 
Tp = 296 MeV: the polarization transfer observables 
have been, in this case, quite instrumental to separate 
the non-spin flip AGDR from the spin-flip SDR in the 
multipole decomposition analysis. The energy differ¬ 


ence between the AGDR and the IAS was determined 
to be Aagdr — Aias = 8.69 ± 0.36 MeV, where the 
uncertainty is claimed to include both statistical and 
systematic contributions. We show this datum by a 
solid (magenta) line in Fig. 4. The other experimental 
measurement has been reported in Ref. [HS ] (Exp2): in 
this case, the ^°®Pb(p, n'yp) ^°^Pb reaction at a beam 
energy of 30 MeV has been used to excite the AGDR 
and to measure its y-decay to the isobaric analog state, 
in coincidence with proton decay of the IAS. The en¬ 
ergy difference Aagdr — Ajas between the AGDR and 
the IAS was determined to be Aagdr — Aias = 8.90 ± 
0.09 MeV. We show this result by means of a dashed 
(blue) line in Fig. 4. Given the error bars, the two work 
provide consistent results. 



FIG. 5: (Color online) The (a) upper and (b) lower pan¬ 
els show the correlations between the neutron-skin thickness 
and either the symmetry energy J at saturation density or 
the corresponding slope parameter L, respectively. The con¬ 
straints provided by the experimental data already shown in 
Fig. 4. 


By comparing the experimental data for the energy 
difference of AGDR and IAS with our correlation line, 
we find that the value of the neutron-skin thickness 
of ^°®Pb is ARpn = 0.212 ± 0.028 fm using Expl, 
and ARpn = 0.195 ± 0.007 fm using Exp2, respec¬ 
tively: these values are indicated by arrows in Fig. 
4. The weighted average of these two results, that 
is, 0.204 ± 0.09 fm for ARpn, is reported in Table I. 
We also compare this result with previous results de¬ 
duced from different experimental methods. Although 





















FIG. 6: (Color online) The values of the slope parameter 
L and symmetry energy J at saturation density extracted 
in the current work are compared with the values extracted 
from other experimental data with several different methods. 


on the high side, our result is compatible with several 
other estimates. In Ref. [sj, [57 1, the authors also 
compared their experimental data with the values of 
T^agdr — T^ias of AGDR and IAS obtained by using 
a fully self-consistent proton-neutron relativistic RPA 
with a family of density-dependent meson-exchange in¬ 
teractions (DD-ME) [^. Including the uncertainty 
both from experimental and theoretical sides, finally 
they found the value of the neutron-skin thickness of 
208pb to be Ai?p„ = 0.216 ± 0.046 fm ± 0.015 fm for 
Expl, and Ai?p„ = 0.190 ± 0.028 fm for Exp2, respec¬ 
tively: these agree well with our results obtained using 
non-relativistic Skyrme energy density functionals. 

We have also extracted the symmetry energy J and 
its slope parameter L at saturation density by using the 
neutron-skin thickness presently obtained. The results 
are shown in Fig. 5. The value for symmetry energy J 
is extracted to be J = 31.8 ± 1.6 MeV (J = 30.9 ± 0.5 
MeV) from Expl (Exp2) at the saturation density, and 
the value for the slope parameter L of symmetry energy 
at saturation density is L = 81.8 ± 17 MeV [L = 71 
± 4 MeV) for Expl (Exp2). The weigthed average of 
these results is J = 31.4 ± 0.5 MeV and L = 76.4 ± 5.4 
MeV. Of course, by making the weighted average one 
reduces the error bars, and this may hide even further 
systematic differences between the experiments and/or 
model dependences. As we mentioned in the abstract, 
the reported errors correspond to a lower-limit estimate 
of the systematic plus experimental uncertainties. 

In Fig. 6, the extracted values of J and L by the 
present analysis are shown together with those obtained 
with other methods. These include: Quantum Monte 
Carlo (QMC) simulations of neutron stars [6^. a nalysis 
of the nuclear binding energies (by FRDM) [^, ener¬ 
gies of isobaric analog states (IAS) 0, proton elastic 
scattering (^°®Pb (p,p)) [1^, pygmy dipole resonances 
(PDR) [^, total dipole polarizability [^, and excita¬ 
tion energy of the isovector giant quadrupole resonance 
[43 |. We should note that the presently extracted value 
of J is consistent with the other values in Fig. 6, with 
small variations. On the other hand, the present value 


of L, although similar to those from the IAS aualysis 
and binding energies from FRDM, is somewhat larger 
than the average value of all other deductions. 

In Ref. [13, the values of L and J were extracted by 
using the same experimental energy difference of AGDR 
and IAS that we have used (Exp2). They have obtained 
J = 32.7 ± 0.6 MeV and L = 49.7 ±4.4 MeV, as shown 
in Fig. 6 with red shaded area. While our result for 
J is consistent with the one obtained in Ref. [13, the 
present (central/ value of L is about 40% larger than 
that of Ref. [13. This may be due to the diffenent 
energy density functionals used iu the present analysis 
and in Ref. [13 , where the RMF Lagrangians of DD-ME 
type were adopted. The inputs for fitting the DD-ME 
and SAMi-J functionals are not exactly the same, and a 
different ansatz for the density dependence is assumed. 
This shows up in, e.g., different values for the nuclear 
incompressibility that turns out to be A'oo = 270 MeV 
in the relativistic case and A"oo=245 MeV for the non- 
relativistic SAMi family. However, a clear explanation 
of this difference is a point that remains for future study. 

We should also notice that in the previous studies 
devoted to the extraction of L and J from giant res¬ 
onances, we have found values of L like 64.8 ± 15.7 
from PDR and 37 ± 18 from IVGQR as reviewed in 
[13 . These values are smaller than the present value: 
the present value is consistent with one of the previ¬ 
ous estimates but not with both of them. On the other 
hand, J is consistent with other extractions from giant 
resonance data. 


V. SUMMARY AND PERSPECTIVE 

In this work, we have studied the correlation of 
the neutron-skin thickness and the energy difference 
Aagdr — Aias of AGDR and IAS in ^°®Pb, by us¬ 
ing a family of effective Skyrme energy density func¬ 
tionals, named SAMi-J (SAMi-m), that are character¬ 
ized by different values of symmetry energy J (effec¬ 
tive mass m*). The calculations have been done within 
a fully self-consistent Skyrme HF plus charge-exchange 
RPA framework. We find a strong linear correlation of 
the energy difference Fagdr — Ajas with the neutron- 
skin thickness Ai?p„ in ^°®Pb. An analytic model has 
been developed to explain the dependence of the ex¬ 
citation energy of AGDR on the neutron-skin thick¬ 
ness, in which it becomes apparent that such excita¬ 
tion energy decreases when the neutron-skin thickness 
increases. We also confirmed that the symmetry energy 
J and the slope parameter L have linear correlations 
with the the neutron-skin thickness within the employed 
Skyrme SAMi-J models. 

Accordingly, we have extracted the neutron-skin 
thickness in ^°®Pb as Ai?p„=0.204±0.009 fm by com¬ 
paring with the corresponding experimental energy dif¬ 
ference of AGDR and IAS. Finally, we have also con¬ 
strained the symmetry energy (J = 31.4 ± 0.5 MeV) 
and its slope parameter {L = 76.4 ± 5.4 MeV) at sat¬ 
uration density by using the value of the neutron skin. 
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Good agreement is obtained in comparing with our new 
results for the neutron-skin thickness and the symme¬ 
try energy J with the values extracted with many dif¬ 
ferent experimental methods. On the other hand, the 
presently extracted L value is somewhat larger than the 
previously obtained values. The reported errors in our 
theoretical analysis correspond to a lower-limit estimate 
of the systematic plus experimental uncertainties. 

The use of the (p, n) reaction to study the AGDR 
can be extended to unstable nuclei due to the progress 
made in the development of new experimental tech- 
niques involving radioactive beams in inverse kinematics 
[73 . . Further experimental efforts on the AGDR in 

other mass regions and/or in long isotopic chains are 
desirable to increase the predictive power of current en¬ 
ergy density functionals and to reduce the model de¬ 
pendence that one deals with when extracting nuclear 


matter properties. This may eventually allow us to bet¬ 
ter constrain the equation of state of asymmetric nuclear 
matter, a landmark for nuclear physics and nuclear as¬ 
trophysics. 
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